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Deterministic lateral displacement (DLD) devices separate micrometer-scale particles in solution based on 
their size using a laminar microfluidic flow in an array of obstacles. We investigate array geometries with 
rational row-shift fractions in DLD devices by use of a simple model including both advection and diffusion. 
Our model predicts novel multi-directional sorting modes that could be experimentally tested in high-throughput 
DLD devices containing obstacles that are much smaller than the separation between obstacles. 
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I. INTRODUCTION 

Deterministic lateral displacement (DLD) is a mechanism 
of particle separation that uses the laminar properties of mi- 
crofluidic flows in a periodic array of posts to sort particles 
based on size. This technique has been shown to differentiate 
between micrometer-sized particles with a resolution in diam- 
eter on the order of 20 nm. The basic sorting mechanism has 
been described for the devices used experimentally: particles 
smaller than a critical radius r c follow streamlines through the 
array while larger particles are systematically 'bumped' later- 
ally during each interaction with a post [1-3]. 

Previous analysis of DLD sorting has focused on predicting 
r c as a function of array parameters, typically the width of the 
gap between posts and the shift of posts between rows. Once 
basic hydrodynamics is included, theoretical calculations of 
r c agree with experimental results within about 5 % [2, 4, 
5]. Inclusion of diffusion in DLD sorting has been described 
using rough estimations [1-3], and in more detailed studies 
that incorporate both microfluidic advection and diffusion to 
calculate r c under a range of experimental conditions [5]. 

Previous analysis of the geometry of the DLD array has 
been limited to the following conventional case. In a given 
row the center-to-center distance between the posts is denoted 
A, see Fig. 1 . The subsequent row of posts is placed at a dis- 
tance aX downstream from the first row. Normally, a is cho- 
sen to be unity, however this is not an essential requirement. 
The posts in this second row are displaced a distance {1/N) X 
along the row, where N traditionally has been an integer. The 
ratio 1 /N is also denoted the row-shift fraction e. In row num- 
ber N + 1 the posts have the same positions as in the first row, 
and consequently the array is cyclic with period 7Y. Due to 
this periodicity of the array and the laminarity of the flow, the 
stream can naturally be divided into TV flow lanes, each carry- 
ing the same amount of fluid flux, and each having a specific 
path through the device [1]. 



For devices with the simple row-shift fraction e = 1/N 
and disregarding particle diffusion, only one critical separa- 
tion size r c is introduced. Spherical particles with a radius 
smaller than r c will move forward along the main flow direc- 
tion through the device, defining the angle 9 = 0. However, 
particles with a radius larger than r c are forced by collisions 
with the posts to move in a skew direction at an angle 9 given 
by tan# = l/(aN). Taking diffusion into account the tran- 
sition from straight to skew motion takes place over a finite 
range of particle sizes [1-3, 5]. 

In this paper we generalize the array geometry by study- 
ing the effects of row-shift fractions different from that of 
the conventional, simple (l/iV)-array. We show in Sec. II 
that by displacing consecutive rows by the rational fraction 
eX = (M/N)X, where M is an integer that is not a divisor of 
N, two new separation modes appear, each associated with a 
distinctive range of particle sizes and separation directions 6. 
Furthermore, to test experimental feasibility of the novel sep- 
aration modes, we introduce in Sec. Ill a model of the DLD 
system reduced to its essential elements: particle trajectories 
interrupted by size-dependent interactions with a periodic ar- 
ray of posts. Utilizing these simplifications, we investigate in 
Sec. Ill A the advection and diffusion of particles in the M/N- 
array geometries, and discuss in Sec. IV possible experimental 
consequences of our novel DLD system. 

In our model of the DLD system described in Sec. Ill we re- 
duce the posts to point-like obstacles in a uniform flow. This 
particular case is currently of interest to researchers looking 
to apply DLD separation to high-throughput microfluidic de- 
vices. Such a reduced post size decreases hydraulic resis- 
tance and thus increases the liquid throughput for a given 
pressure difference applied along the device. One promising 
method to create such devices is to use arrays of semiconduc- 
tor nanowires [6] in a microfluidic channel. 



II. BASIC THEORETICAL ANALYSIS 



'Current address: Biomedicial Engineering Division, Oregon Health and Sci- 
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The introduction of a non-simple row-shift fraction e = 
M/N in the DLD system is first discussed in Sec. II A for 
the specific case of M/N = 3/8, since all the novel separa- 
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FIG. 1: An array of posts (marked by black dots) with period N = 8 
and a shift of M = 3 flow lanes per row, i.e., a row-shift fraction 
e = 3/8 = 0.375. The flow v is directed along the x-axis from 
top to bottom. The dashed lines indicate the four possible separation 
directions. First, the two well-known modes: the straight mode A 
with 6 a = 0° for particles of radius r with r < (1/N)X (small 
dark gray circles starting in flow lane 7? = 1), and the maximal 
displacement mode C with 9 C = 20° for (3/N)X < r (large light 
gray circles starting in flow lane R = 4). Additionally, the two novel 
separation modes: one B+ towards the right with angle 8b + = 2.4° 
for (1/N)X < r < (2/N)X (small light gray circles starting in flow 
lane R = 2), and another B- towards the left with angle 0b_ = 
-7.1° for {2/N)X < r < {3/N)X (large dark gray circles starting 
in flow lane R = 3). The solid vertical lines indicate the flow lanes 
of width X/N, while ai and a2 are the lattice vectors, and a is an 
aspect ratio. 



tion modes are present in that device geometry. Fig. 1 shows 
the principle of the fractionally displaced DLD array leading 
to multi-directional separation of particles of different sizes. 



As was the case for the simple row shift fraction 1/N, the ra- 
tional row-shift fraction M/N also naturally leads to N flow 
lanes, each carrying the same amount of fluid flux. In this sec- 
tion, all particles are assumed to follow these flow lanes unless 
bumped by an interaction with a post. However, in contrast to 
the traditional DLD geometries, now the posts are displaced 
M flow lanes instead of just a single flow lane when passing 
from one row of posts to the next. 

In Sec. II B we analyze this more general case of M/N- 
arrays, where the integer row-shift M and the integer array- 
period N have no common divisors. 

A. The specific row-shift fraction 3/8 

First we consider the explicit choice of parameters given 
in Fig. 1, namely, a period N = 8, and a row-shift of 
M = 3 lanes in the y-direction e y , i.e., a row-shift fraction 
of e = 3/8. The flow v is in the x-direction e x . For simplic- 
ity, we employ the most simple model where all flow lanes are 
assumed to have the same width X/N, and where the particles 
are not subject to Brownian motion. The analysis can straight- 
forwardly be extended to take the different widths of the flow 
lanes [2] as well as diffusion [5] into account. 

The analysis is most easily carried out by considering 
spherical particles of increasing radius r. As the rows in Fig. 1 
are shifted to the right, it is natural to choose the starting point 
of a given particle to be directly to the right of a post, placing 
the particle's center in flow lane R = 1, 2, 3, or 4 according to 
size. 

For the smallest particles with r < X/N, labeled A in 
Fig. 1, we obtain a path corresponding to the familiar so- 
called zigzag path defined in Ref. [1]. Due to the point- 
like nature of our obstacles, the path is a straight line, indi- 
cated by the dashed vertical line in Fig. 1. The path angle is 
9 A = arctanO = 0°. 

For the next set of particles with X/N < r < 2X/N, (B + in 
Fig. 1), we note that they are not affected significantly by pass- 
ing the second rows of posts. The displacement of (M/N)X 
is larger than the size of the particle. By simple inspection 
we find that the particles interact with a post in the fourth row 
leading to a bump of one lane width to the right. This bumping 
brings the particles back to a position just right of a post, and 
we have identified a new separation mode, B + . The direction 
of mode B + can be characterized by the integers 

p = the number of rows after which the bumping (la) 
pattern repeats itself, and 

q = the number of flow lanes that the particles are (lb) 
bumped to the right while traveling through p 
rows. 

Here, with p = 3 and q — 1 and the array parameters indicated 
in Fig. 1, the path angle of mode B + is found to be 0b + = 
arctan [l/(a x 3 x 8)] = 2.4°. Here and in the following we 
choose the aspect ratio a = 1. 

For the third set of particles with 2X/N < r < 3X/N, 
marked as £?_ in Fig. 1, we note that they collide with a post in 
the second row and are bumped two lanes to the left. After two 
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TABLE I: List of separation radii r 


and angles 6 


as a function of the 


integer 


array parameters N, M, p and q for a = 


1. 


AT 71//" 


mode 


particle radius in 


q/p 


separation angle 






units of lane width 


u — arctan [^^J 


5 2 


^4 


< r < 1 


0/1 


u.u 




_R_ 


1< r < 2 


-1/2 


-5.7° 




(7 


2 < r < 2.5 


2/1 


21.8° 


7 2 


^4 


< r < 1 


0/1 


0.0° 




R 


1 < r < 2 


-1/3 


-2.7° 




c 


2 < r < 3.5 


2/1 


15.9° 


7 3 


A 


< r < 1 


0/1 


0.0° 




B- 


1 < r < 3 


-1/2 


—4.1° 




C 


3 < r < 3.5 


3/1 


23.2° 


& ^ 

o o 


A 
I\ 


< r < 1 


0/1 


0.0° 




R , 


1< r < 2 


1/3 


9 A° 
Z.I 




B_ 


2 < r < 3 


-2/2 


-7.1° 




<j 


3 < r < 4.0 


3/1 


20.6° 


9 2 




< r < 1 


0/1 


0.0° 




B- 


1 < r < 2 


-1/4 


-1.6° 






2 < r < 4.5 


2/1 


12.5° 


9 4 


A 


< r < 1 


0/1 


0.0° 




B- 


1 < r < 4 


-1/2 


-3.2° 




c 


4 < r < 4.5 


4/1 


24.0° 


10 3 




< r < 1 


0/1 


u.u 




R 




1 l"\ 


-1.9° 




c 


3 < r < 5.0 


3/1 


16.7° 


11 2 


A 


U / <s. 1 


n/i 

U/ 1 


0.0° 




B- 


1< r < 2 


-1/5 


-1.0° 




C 


2 < r < 5.5 


2/1 


10.3° 


11 3 


A 


U V v <, 1 


n /i 

U/l 


0.0° 




B+ 


1 < r < 2 


1 /4 


1.3° 




B- 






-3.5° 




C 


3 < r < 5.5 


3/1 


15.3° 


11 4 


A 


< r < 1 


0/1 


0.0° 




a + 


1 < r < 3 


1/3 


1.7° 






3 < r < 4 


-3/2 


-7.8° 




c 


4 < r < 5.5 


4/1 


20.0° 


11 5 


A 


< r < 1 


0/1 


0.0° 




B- 


1 < r < 5 


-1/2 


—2.6° 




C 


5 < r < 5.5 


5/1 


24.4° 


12 5 


A 


< r < 1 


0/1 


0.0° 




B+ 


1 < r < 2 


1/5 


1.0° 




B- 


2 < r < 5 


-2/2 


-4.8° 




C 


5 < r < 6.0 


5/1 


22.6° 



rows, the particles are again bumped two lanes to the left, and 
we have identified another new separation mode, £>_ . Given 
this period p = 2 bumping of q = — 2 flow lanes (where minus 
indicates displacement to the left), the path angle of mode _B_ 
is found to be 8 B _ = arctan [ - 2/ (a x 2 x 8)] = -7.1°. 

Finally, the fourth set of particles (with 3 A/TV < r) is con- 
sidered, shown as the large light gray circle in Fig. 1. Since 
3 A/TV equals the row-shift eA, these large particles collide 
with a post in each row (p = 1) where they are bumped 
q = M = 3 lanes to the right. This is the conventional max- 
imal displacement mode C [1]. As a result the path angle for 
mode C here is found to be 9c = arctan [3/ (ax 8)] = 20.6°. 



B. General row-shift fractions M/N 

In the general case of a DLD device with period N and a 
row-shift of M flow lanes, it is useful to introduce the floor 
function [x\ of x, which gives the largest integer smaller than 
or equal to x, eg., [8/3J =2 and [10 /3J = 3, and the ceiling 
function \x~\ of x which gives the smallest integer larger than 
or equal to x (see also the definitions given at Ref. [7]). 

Using the notation in Fig. 1 , the flow lane R occupied by the 
center of the particles can be expressed in terms of the particle 
radius r as R = \rN/ A] , so that R = 1,2,3, ... , \N/2] 
forO<r<A/2. 

Two cases are straightforward to analyze. For small radii 
with R = 1, the particles will follow the streamlines without 
any systematic net lateral displacement, i.e., a mode A in the 
direction given by 



ae. 



and forming the path angle 9a with the x-axis, 
6^ = 0, R=l. 



(2) 



(3) 



For large radii with M < R < |"^] , the particles collide with 
the posts and are bumped M flow lanes to the right in each 
row, but they do not get stuck between the posts; this is mode 
C. The path is directed along the direction tc given by 



tc 



M 



and forming the path angle 9q with the a>axis, 



9c = arctan 



M 



aN 



M <R< 



TV 
~2 



(4) 



(5) 



In a given M/N array, modes with larger sorting angles are 
excluded because of the post spacing in the y-direction: parti- 
cles with radius r > A/2 are unable to fit between the posts. 

If the particles are small enough to pass the second row 
without getting bumped to the right, but too large for mode A, 
1 < R < M, their trajectories fall into one or two B modes. 

As a particle is convected through the array, a post will ap- 
proach the particle from the left in steps of M flow-lanes per 
row the particle advances, hence the use of modulus M arith- 
metic in the following analysis. 

If (JV mod M) < R < M the particle will hit the post 
with its center to the left of this obstacle and will therefore 
enter mode £>_ where it is displaced to the left with a period 
p_ = . This is most readily seen by starting the analysis 
with a particle position just left of a post. A particle with 
(N mod M) < R < [{N - R + 1) mod M] will bump 
left after p_ = [-p-J rows and will again be in a position 
just left of a post. The small particle in mode £?_ of Fig. 2 
is an example of this behavior. Slightly larger particles with 
[(N - R + 1) mod M] < R < M will bump right after 
p = [ N ~m +1 \ rows. Since we are only considering particles 
with R < M, this displacement will always be less than M 
flow lanes, and the particle is therefore bound to bump left on 
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FIG. 2: An array with period N — 10 and a shift of M = 3 flow 
lanes per row, giving a row-shift fraction e = 3/10 = 0.3. The 
flow v is directed along the x-axis from top to bottom. Here there 
are three sorting modes, delimited by two critical radii. Mode A for 
particles of radius r with r < r c i = (1/10)A (shown on the far left), 
and mode C, the maximal displacement mode for (3/10)A = r C 2 < 
r < (1/2) A (large light gray circles). A novel mode B- displaces 
particles with r c i < r < r C 2 towards the left (large dark gray and 
intermediate, light gray circles). 



the post in the following row, i.e., after a total of p_ = [^J 
rows. The large _B_ mode particle in Fig. 2 is an example of 
this behavior. 

The trajectories in mode £?_ have period p_ = . The 
number q_ of lanes bumped after passing these p_ rows is 
<7_ = Mp- — N < 0. The path is directed along the direction 
ts_ given by 



P- 

1V1 

q- = Mp- - N < 0, 
forming the path angle 0g_ with the x-axis, 



ctp- e x 

N ' 
M_ 
Mp- 



9b_ = arctan 



ap-N 



< 0. 



(6a) 

(6b) 
(6c) 

(7) 



If 1 < R < (N mod M) the particle will enter mode B + 
where it is displaced to the right with a period p+ = [ N ^ R \ ■ 
To realize this it is natural to start the analysis with the particle 
just right of a post. Again, a post will approach the particle 
from the left in steps of M lanes as the particle moves through 
the array. A particle with 1 < R < [~4pl will follow the 
flow forp + = [ N ^j R \ rows and then bump right. If [4^1 < 
R < (N mod M) the particle will bump left already in the 
second row of posts. The particle is now in a position just 
left of a post. However, since it is not large enough to follow 
the B- path, it will bump right when it meets the post after 
P+ = L^J rows. 

The trajectories in mode B+ have period p + = \_ N ^ R \ ■ 
After p + rows the particles will get bumped q + flow lanes to 
the right given by q + = Mp + — N > 0. The path is directed 
along the direction t b + given by 



t B+ = ap+ e K + — e. y , 



P+ 



N + R 
M 



q+ = Mp+ -N>0, 
forming the path angle 9b + with the .T-axis, 



arctan 



ap+N 



> 0. 



(8a) 

(8b) 
(8c) 

(9) 



In terms of the flow lane number R, the criteria for the four 
different displacement modes can be summarized as follows 



mode A, 
mode B + , 
mode £>_, 
mode C, 



ifi?=l (10a) 

if K R < (N mod M) (10b) 

if (N mod M) <R<M (10c) 

if M < R< [f ]. (lOd) 



Note that mode B + vanishes if (N mod M) = 1, 
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in. MODEL AND IMPLEMENTATION 

The following model is established to numerically test the 
sorting behavior of a particular M/N DLD array and take into 
account the effect of particle diffusion on sorting behavior, as 
discussed below. We treat the device as a periodic array of 
zero-radius posts with the geometry shown in Fig. 2. This 
N = 10, M = 3 geometry, with a row-shift fraction given by 
e = 3/10, exhibits the three modes shown in Table I, includ- 
ing a novel sorting mode, _B_ . 

We assume the array to be infinitely deep so that the flow 
field is two-dimensional and independent of the z-direction. 
Consistent with the infinitesimal size of the posts, the liquid 
flow through the device is assumed to be uniform with ve- 
locity v = v e x along the ir-axis. Thus, our model does not 
describe Taylor-Aris dispersion, which in real systems with 
finite-sized posts would be induced along the x-direction by 
a combination of transverse diffusion and transverse veloc- 
ity gradients [8]. The particles only interact with the posts 
through a hard-wall repulsion and any effect of the particles 
on fluid flow is neglected. The particle-post interaction ex- 
cludes the center of a particle with radius r from a circular 
region of the same radius around the point-sized post. In ad- 
dition to being moved by the fluid and interacting with the 
posts, each particle has a diffusion coefficient D given by the 
Einstein relation 

D(r) = ^, (11) 

where k^T is the thermal energy and r\ is the viscosity of the 
water. For the calculations below we have chosen the fol- 
lowing experimentally relevant parameters: For water at room 
temperature k B T « 4 x 10~ 21 J and r\ w 10~ 3 Pa s, and for 
the geometry the post separation is A = 10 |im and particle 
radii in the range 0.5 |im < r < 4 p.m. 

A final basic assumption of our model is that all time de- 
pendence in our model is implicitly given by the advective 
flow speed v. For particles starting at the entrance of the de- 
vice at x = the time t is given through its x-coordinate 
as t = x/v. The model therefore allows all the relevant dy- 
namics of an ensemble of many particles to be described by a 
continuous concentration distribution c(x 1 y) with some given 
initial distribution c(0, y) at the entrance of the DLD device. 
Given c(0, y) the time-evolution of the distribution consists of 
calculating c(Ax, y) after convection to x = Ax. By follow- 
ing the evolution of c(x, y) as the distribution interacts with 
posts and responds to thermal forces, our model can identify 
the basic modes of transport in an array of posts and the effect 
of diffusion on this transport. 

The initial distribution c(0, y) is given by a box distribution 
of width A (although a narrow distribution is used in Fig. 3 for 
visual clarity), and the distribution c(Ax, y) is calculated from 
the previous distribution c(0, y) taking into account its inter- 
actions with the posts as well as the diffusion equation. The 
entire distribution c(x, y) is evaluated by iterating the follow- 
ing procedure: 

1 . Upon encountering a row of posts, the distribution for 
particles of radius r is set to zero in regions with a dis- 



tance smaller than r to any post, and the corresponding 
number of particles is then added to the distribution in 
the adjacent pixels to maintain the total number of par- 
ticles (see Fig. 3). 

2. The distribution c(x, y) is subsequently evolved in ac- 
cordance with the diffusion equation, with the diffusion 
coefficient given by Eq. (11), 



dc _ d 2 c 
dx dy 2 ' 



(12) 



employing the implicit time t — x/v set by convec- 
tion along the x-direction, and using the Fourier co- 
sine transformation in the transverse y-direction as de- 
scribed below. 

The computation uses a finite array of width w = 10A, i.e. 
containing 10 posts, and the row separation is again taken to 
be equal to the post separation, i.e. a = 1. The array with 
width w is discretized in y into n max = 10 4 pixels of size 
Aw x Ax with Ax = Aw = w/n max . 

The discrete Fourier cosine transformation C(x, k n ) of the 
distribution c(x, y) then takes the form 



C(x, k n ) 



2 — 5o,n 



w 



c(x,y) cos(k n y) dy, (13) 



where k n is given by 

k n = n, 
w 



71 — 0,1,2,..., VT-niax — 



W 

Aw' 



(14) 



By direct inspection we find the well-known result from the y- 
dependent diffusion equation in Fourier space that, during the 
time step At = Ax/v, C(x, k n ) evolves into C(x + Ax, k n ) 
as 



C{x + Ax, k n ) = C(x, k n ) exp 



-Dkl 



Ax 



(15) 



By the inverse Fourier cosine transform we can therefore 
write the distribution at row x + Ax in terms of that at row 

x as 



c(x + Ax,y) = ^ C(x,k n ) cxp 



n=0 



Dkl— cos(k n y). 

(16) 

which by construction automatically respects the boundary 
condition that no particles can diffuse beyond the edges of 
the array. The evolution of the distribution due to diffusion is 
computed at each row of pixels after the effects of posts on the 
distribution have been taken into account. 

To elucidate the sorting mechanism in the absence of ther- 
mal forces, calculations were also done with diffusion coeffi- 
cient D = 0, in which case c(x, y) evolves only according to 
the interaction of the particles with the posts. Results of these 
calculations are shown as D = in Fig. 4. 

While c(x, y) is the calculated distribution at a given time 
and position in the array, the set of all c(x, y) also represents 
the steady-state distribution of a stream of particles entering 
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FIG. 3: [Color online] Composite image of numerically calculated 
spatial distributions for a device with N = 10 and M = 3. (a) 
Results for D — for particles with radii r < r c i (mode A), r c i < 
r < r c2 (mode B) and r > r C 2 (mode C). (b) Particles, with the 
same radii as in (a), moving through the array with a flow speed of 
v — 100 um/s, including the effect of diffusion. Broadening of all 
distributions due to diffusion can be seen and particles in mode C are 
sorted less efficiently. Initial spatial distributions here are the same 
for all particle radii, and a narrow initial distribution is used for visual 
clarity. 



an array of obstacles and moving constantly through the array, 
as seen in Fig. 3. 

The calculations were done using Matlab on a personal 
computer and a 64-bit dual processor workstation. 



A. Results 

1. Three transport modes in the 3/10-array 

The existence of the novel sorting mode B-, as well as 
the two modes A and C previously described in DLD liter- 
ature are confirmed by applying our numerical model to a 
range of particle sizes advected through the 3/10-array. As 
the particle distributions move through the array, their trajec- 
tories form three modes A, B and C, according to two critical 
radii, r c \ and r C 2, see Fig. 3(a). Our calculations reproduce 
the two known modes: the 'zigzag mode' A, in which there 
is no average displacement from the direction of flow, and 
the 'bumped mode' C, in which particles are bumped later- 
ally in every row. These two modes are most clearly seen in 
Fig. 3(a), where the distributions are calculated without diffu- 
sion. In mode A, where r < r c i = (1/N)X, particles may 



interact with the posts, but no net lateral displacement is ac- 
complished. Mode C is characterized by a displacement equal 
to the shift [M/N)X for every row the particles pass through. 
In the novel mode particles of size r c \ < r < r C 2 interact 
with posts more frequently than in mode A but less frequently 
than in Mode C, as described in Sees. II B. The 3/10 array 
used here clearly exhibits the lone £>_ mode shown in Table I 
for these array parameters. It is important to note that mode 
£>_ vanishes in the conventional case M = 1, and all particles 
smaller than the critical radius r C 2 move along the direction of 
flow. 

The directions and tc of the conventional modes A and 
C are given directly by Eqs. (2) and (4) for a = 1: 



U = e 2 



3 

10 ^ 



(17a) 
(17b) 



while the direction t b_ of mode B_ is found through the path 

period v- = \ E ir\ = = 3 and the lane shift 1- = 

p_M - N = 3 X 3 - 10 = -1, and thus 



t B _ =3e x - — e v . (17c) 
The corresponding displacement angles become 

0a= o.o°, oA <r< iA ; ( i 8a) 
e B _ = -i.9°, 1 ^< r < 3 ^ ( 18b ) 

C = 16.7°, 3A< r<5 A. (i 8c) 

The array parameters used here can be translated into those 
used in DLD literature [1-3], simply by setting M = 1. 



2. Effect of diffusion on sorting 

The effect of diffusion on the sorting of particles is shown 
in Fig. 4. The angles shown are measured between v and the 
lateral displacement of the center of mass of the distribution 
for each particle size after 10 rows of posts for high and low 
flow speeds. We can estimate speeds at which diffusion be- 
comes negligible by comparing the time it takes a particle to 
be advected along the x-direction from one row to the next, 
X/v, to the time it takes a particle to diffuse transversely in 
the y-direction to reach a position where it would be bumped, 
2D/(r — r C 2) 2 . For high flow speeds, 



v > 



2DX 



(r - r c2 ) 2 ' 



(19) 



diffusion can be neglected, and the transitions between the 
sorting modes are sharp, as seen in the D = case. Note 
that this velocity diverges as the particle size approaches the 
critical radius r c i\ in this limit the displacement needed for a 
particle to change sorting directions goes to zero. Within the 
spatial resolution of this work (1 pixel =10 nm), the particles 
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FIG. 4: Sorting angles 9 calculated for a device with N = 10 and 
M = 3 as described in Sec. Ill with r c i = 1 um and r C 2 = 3 um. The 
initial position corresponds to the center of mass of the initial distri- 
bution between two posts and the final position corresponds to the 
center of mass after 10 rows. 9 is plotted here versus particle radius 
r with (open red circles) and without (filled blue circles) diffusion 
with flow velocity v = 10 um/s. The negative sorting angles for 
r c i < r < r C 2 indicate the presence of mode B_ for this array. Inset 
shows the sorting angle 9 around r — r C 2 = 3 um for a range of flow 
velocities (same y-axis range). Diffusion blurs the sharp transition 
between the sorting modes, as discussed in Sec. IV B. 



closest in size to the critical radius will still be sensitive to 
diffusion at flow velocities below 10 mm/s. As flow speeds 
decrease, particles have more time to diffuse transversely as 
they move through the array, and the effects of thermal mo- 
tion on sorting are seen more clearly. Transverse diffusion of 
particles along the y-direction tends to shift the center of mass 
of the distribution c(x, y) towards the midpoint between posts. 
This means that particles with r < r C2 = |^ are more likely 
to be shifted to higher sorting angles. However, in the regions 
between rows, diffusion allows particles to move transversely 
away from the path that would normally be 'bumped' by a 
post, decreasing their sorting angle. These two effects of dif- 
fusion are responsible for the smoothing of the angle versus 
radius curves for slower flow speeds in Fig. 4. The calculated 
values for 8 b - are in good agreement with the value predicted 
in Eq. (18a), but for particles with r > r C2 = |{y, the finite 
width of the initial distribution and the relatively short array 
size (10 rows) reduce the calculated values for 9c from the 
predicted value by about 15 %. The small variation in sorting 
angle with radius for modes £>_ and C for D = in Fig. 4 is 
mainly the result of the two end points used to define the angle 
being not exactly equivalent: the position of the second, but 
not the first end point varies continuously with bead size, and 
so the presented angle varies with bead size. Secondly, since 
the number of rows is not divisible by the periodicity of mode 
£?_, an additional small error is introduced. These deviations 
should vanish for simulations with larger numbers of rows. 



IV. DISCUSSION 

A. The novel sorting mode and its relation to kinetically 
locked-in transport 

DLD devices have thus far been made with a fixed flow 
direction and almost exclusively with M = 1. However, in 
theoretical work studying transport through periodic potential 
landscapes, the direction of the applied force is varied for a 
fixed array geometry and the transport direction is calculated 
[9-12]. To calculate the correspondence between varying the 
array parameters M and N used here and changing the flow 
direction in a fixed array as in [9-12] is cumbersome, but for 
a range of flow directions near tg_ = 3ea; — -X e y , the angles 
to the flow direction 9b_ and 9c vary as the flow directions 
change, but the relative angle between them, 9c — 8b_, re- 
mains a constant defined by the array. The angle between 
modes £?_ and C is insensitive to small changes in flow di- 
rection for v near (in this case) 3ea; — e y . 

This insensitivity to flow direction is an example of a 
plateau in a so-called 'devil's staircase' : transport through a 2- 
D periodic potential is independent of the flow direction near 
small integer lattice vectors [9]. In this case the lattice vectors 

e„, and the two 



J_ e 

to c v 



= 3ai -a 2 



are ai = tc — e x + (3/10) e v and a 2 
close-lying flow directions are ts_ = 3e T — 
and ai = t c . 

The interplay between lattice directions and applied forces 
has been documented extensively in the literature of kineti- 
cally and statistically locked-in transport. Of interest in the 
present context is that many numerical simulations of trajecto- 
ries through various two-dimensional periodic potentials have 
been done to study these and other phenomena, including sort- 
ing of particles [9-12]. 

The interaction between posts and particles that we have 
chosen simplifies DLD to a ID distribution that evolves in 
time. This allows the effects of diffusion to be easily incorpo- 
rated into our modeling of the dynamics of the distribution of 
particles. Also, the particular interaction between point-sized 
posts and finite-sized particles depends only on particle size, 
an analysis that seems to be absent from the literature. 



B. Diffusion, detectability and experimental possibilities 

A clear difference between the results in Fig. 4, based on 
zero-sized posts, and those reported in the literature, based 
on finite-sized posts, is that the critical radius (defined as the 
inflection point of the angle vs. radius graph near r = r c ), 
decreases for lower flow velocities in Ref. [1], whereas our 
simulations show a critical radius that is essentially constant. 
When particles have more time to diffuse laterally in reported 
experimental data, ones that previously followed the 'zigzag' 
path follow something closer to Mode C but not the other way 
around. We have identified the difference in size of the posts 
as the primary basis for the difference in symmetry. In the 
gap between the posts, only beads smaller than r c can change 
modes (from A to C) whereas beads larger than r c cannot 
change modes because of steric hindrance. Diffusion between 
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posts is thus asymmetric. On the other hand, between rows 
all beads can change modes equally well so that the effect 
of diffusion is symmetric. This result is most clearly seen in 
two cases: (f) with sufficiently large posts and small spacing 
between rows, diffusion between posts dominates leading to 
asymmetry, and (if) with our needle-like posts, instead diffu- 
sion between the rows dominates, leading to symmetry be- 
tween small and large particles. In devices with large round 
posts such as those in Ref. [1], the flow streams are narrower 
in the gap between the posts than in the region between the 
rows making the asymmetric diffusion even more pronounced. 
The symmetry about r C 2 shows that sorting in this model is ro- 
bust against changes in flow velocity. 

As discussed in Sec. IV A, there is no difference between 
modes A and B- when the flow is directed along the lattice 
direction 3a i — a2, which is equivalent to a conventional array 
with M = 1 and A = 3, instead of along e x . Also, while 
mode S_ for the 3/10 array shown in Fig. 2 is directed away 
from mode C, the mode B + discussed in Sec. II A is deflected 
away from v towards mode C. The absence of modes _B_ and 
B + in previous analyses of DLD experiments stems from the 
use of tilted square arrays with flows chosen such that M = 1 
or more general arrays that are still limited to simple row- 
shifts l/N. In these cases, Modes A and B are the same: 
they both go along the direction of flow. Interestingly, in their 
paper Ref. [2], Inglis et al. mention that they are studying 
simple row-shift fractions e = l/N, with N being an integer, 
but they do not comment on the data points in their Fig. 2 that 
clearly have 6^ l/N. 

Experimental detection of mode B requires that the distri- 
butions of modes A and B must be spatially separated. The 
numerically calculated distributions shown in Fig. 5 exhibit 
four qualitative regimes that could be observed in an experi- 
ment to detect the presence of particle transport in mode B. 

(a) . At very high flow speeds, corresponding to D = 

in the numerical data, the three modes are completely 
separated because each distribution is very narrow. In 
this regime, arbitrary spatial separation can be achieved 
simply by running the particles through a longer array. 

(b) . At high intermediate flow speeds, the distributions have 

widened due to diffusion, but modes A and B are 
clearly distinguishable, despite some overlap. 

(c) . At low intermediate flow speeds, modes A and B over- 

lap enough to prevent resolution of two separate dis- 
tributions. This regime is relevant to DLD device de- 
sign because it would be experimentally observed as an 
anomalous, asymmetric broadening of the distribution 
associated with the 'zigzag' path. 

(d) . At low flow speeds, distributions from modes A and B 

are completely overlapping and it may even be difficult 
to differentiate them from mode C. 

Experimental realization of the regime investigated in this 
model would require arrays made with very small posts to 
minimize hydrodynamic effects on particle trajectories. This 
also corresponds to a reduction in hydrodynamic drag, which 
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FIG. 5: [Color online] Distributions of three particle sizes: r — 0.90 
(green dashed line), 2.00 (blue), and 3.08 p.m (red) after transport 
through ten rows of the N = 10, M — 3 array. The total number 
of particles is the same in each case and each initial distribution (not 
shown) is a square distribution with a narrow width centered on y = 
30 |im. (a) No diffusion, (b) With diffusion and v = 1000 p.m/s. 
(c) With diffusion and v — 100 um/s. (d) With diffusion and v = 
10 (im/s. Panel (b) shows a case where mode B_ could be detected 
experimentally. For the lower speed in panel (c), modes A and B_ 
cannot be resolved, but the combined distribution is broader than 
mode A alone. For the even lower speed in panel (d) the distributions 
of particles in modes A and £>_ are each wider than the separation 
between them and the two modes are completely unresolvable. 



is beneficial for researchers seeking to increase fluid through- 
put of devices. 

As can be seen in Fig. 3, the angle 0b is small compared 
to 9c ■ In order to differentiate between particles traveling 
in modes A and B, size dispersion of beads must be consid- 
ered in addition to broadening due to diffusion. Commercially 
available polystyrene beads used in DLD experiments typi- 
cally have size distributions with widths of less than ±10 %. 
This then requires choosing particles whose size distributions 
are separated by more than 10 %, such as those shown in 
Fig. 5, or the use of a DLD array to create a sufficiently nar- 
row size distribution. If hydrodynamic effects or limitations 
on flow velocity in a particular experiment prevent the novel 
sorting mode from being completely resolved, it may still ap- 
pear as an asymmetric broadening of the distribution of seem- 
ingly undeflected particles, as in Fig. 5. 

In general, the separation angles for a given M/N-array can 
be made larger to the extent that the aspect ratio a can be 
made smaller without risking clogging of the largest particles. 
By consulting Table I, it can be seen that the novel separation 
angle of the 3/10-array is one of the smaller B angles, and 
also, the 3/8-, 3/1 1-, 4/1 1- and 5/12-arrays offer both the B + 
and the £>_ modes. 
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V. CONCLUSIONS 

We have identified novel sorting modes in a model of trans- 
port through a DLD device characterized by row-shift frac- 
tions M /N. Our simple model also reproduces key features of 
DLD arrays, including sorting based on size and the blurring 
of cutoffs between modes due to diffusion. Even if not com- 
pletely resolved, the novel sorting mode has the potential to 
increase spatial broadening of 'zigzag' particle distributions. 
In order to avoid this broadening, adjustable DLD arrays could 
use variable spacing while maintaining a fixed M = 1 geom- 
etry, such as in Ref. [4], or tune flow angles to exactly repro- 
duce the M = 1 condition across a fixed obstacle array using 
techniques such as in Ref. [13]. Our simulations indicate that 
using needle-like posts decreases the shift in critical size due 



to diffusion that has been observed in devices where the post 
separation is on the same scale as the post diameter. Further- 
more, the use of more general array geometries and simplified 
fluid dynamics links this work to the field of kinetically locked 
transport phenomena. 
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